Sum-rules and bath-parametrization for quantum cluster theories 
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We analyze cellular dynamical mean-field theory (CDMFT) and the dynamical cluster approx- 
imation (DC A). We derive exact sum-rules for the hybridization functions and give examples for 
DMFT, CDMFT, and DCA. For impurity solvers based on a Hamiltonian, these sum-rules can be 
used to monitor convergence of the bath-parametrization. We further discuss how the symmetry 
of the cluster naturally leads to a decomposition of the bath Creen matrix into irreducible compo- 
nents, which can be parametrized independently, and give an explicit recipe for finding the optimal 
bath-parametrization. As a benchmark we revisit the one-dimensional Hubbard model. We care- 
fully analyze the evolution of the density as a function of chemical potential and find that, close 
to the Mott transition, convergence with cluster size is unexpectedly slow. In two dimensions we 
find, that we need so many bath-sites to obtain a reliable parametrization that Lanczos calculations 
are hardly feasible with current computers. For such large baths our symmetry-adapted approach 
should prove crucial for finding a reliable bath-parametrization. 

PACS numbers: 71.10.-w,71.27.+a,71.10.Fd,71.30.-|-h 



I. INTRODUCTION 

Strongly correlated materials are characterized by 
the interplay of kinetic energy and sizable short-range 
electronic repulsion, which cannot be described with 
single-particle approaches or standard perturbative the- 
ories. The Dynamical Mean Field Thcoryii^ (DMFT) 
has proven extremely powerful, in particular for the de- 
scription of the correlation-driven Mott transition. Yet 
DMFT is strictly local. A number of recent devel- 
opments are aimed at overcoming this limitation and 
are greatly contributing to the understanding of the 
physics of strongly correlated systems i^i^i^iii^i^ In quan- 
tum cluster theories^ the /c-dependence is introduced by 
considering a small number of sites, instead of the sin- 
gle correlated site of DMFT, and embedding them in a 
bath, i.e., a dynamical mean- field host determined sclf- 
consistently. This is a good approximation when the 
self-energy is reasonably localized, which can be sys- 
tematically improved by considering larger and larger 
clusters.— 

Cluster extensions of DMFT are not unique. Here 
we consider the two main flavors, cellular dynamical 
mean-field theory^ (CDMFT) and the dynamical cluster 
approximation^ (DCA). Finite temperature Monte Carlo 
provides an efficient cluster solver, however, to obtain 
spectra on the real-frequency axis, data for imaginary 
time must be analytically continued, usually applying 
maximum entropy. If the spectrum has structures on 
a small energy scale, this approach can lead to problems. 
It is then natural to explore Hamiltonian-based solvers, 
like the Lanczos method, which directly give result on 
the real- frequency axis. 

Lanczos is extensively applied to DMFT and results 
are reliable and very accurate. The critical step in such 
calculations is the fitting of the bath dcgrces-of-frcedom. 
In cluster methods this step gets more involved. This 



is one of the reasons why it found few application for 
cluster methods: To the best of our knowledge, it has 
so far only been used for CDMFT calculations of very 
small clusters. Our aim is therefore to give a systematic 
formulation, exploiting symmetries to arrive at an opti- 
mal parametrization of the bath. In addition we derive 
exact sum-rules for the bath Green's functions. As an 
application, we investigate the one-dimensional Hubbard 
model with CDMFT in detail,— i^i^^ paying particular 
attention to the convergence with the size of both bath 
and cluster. 

The paper is organized as follows: To fix the notation, 
in section |TT1 we give a unified formulation of CDMFT 
and DCA in a formalism using a Hamiltonian solver. 
In section IIIII we derive general sum-rules for the hy- 
bridizations and give examples for DMFT, CDMFT, and 
DCA. In section HV] we discuss how the symmetry of the 
cluster naturally leads to a decomposition of the bath 
Green matrix into irreducible components, which can be 
parametrized independently. We give an explicit recipe 
for finding the optimal bath-parametrization and discuss 
how this approach relates to the technique of cluster 
replica. In section |V] we analyze how the cluster ap- 
proaches work for the Hubbard chain and carefully reex- 
amine the evolution of the density as a function of chem- 
ical potential. In section IVll we give our conclusions and 
an outlook. 



II. METHOD AND NOTATION 

We consider the Hubbard model 



(1) 



To fix the notation we briefly sketch the self-consistency 
loop for cellular DMFT (CDMFT) and the dynamical 



cluster approximation (DCA) using, e.g., exact diago- 
nalization as impurity solver. Let N^ be the number of 
cluster-sites, N^ the number of bath-sites. For simplicity 
we suppress spin-indices. 

Given an Nc x Nc bath Green matrix Q , 

1. fit parameters of an Anderson model with Nh bath- 
sites 



g-p}^^{uj) ^u + ^l-n,-T[u- E]-ir 



irt 



(2) 



to ^^ , where T is the Nc x iV;,-dimensional hy- 
bridization matrix, and E the N}, x A'^^-dimcnsional 
bath-matrix. He is specified below, 

2. solve the A^c + -^&-site Anderson model i^And (spec- 
ified below) to obtain the N^ x N^ cluster Green 
matrix G^, 

3. get the cluster self-energy matrix 

s,H = g-iH-G-iM, (3) 

4. calculate the local Green matrix for the cluster by 
integrating over the reduced Brillouin-zone of the 
cluster 



G{lo) = dkiuj + ^i- H(k) - Sc(cj) 



(4) 



where H(k) is the single-electron part of the of the 
Hubbard Hamiltonian ^ in the reduced Brillouin- 
zone of the cluster, 

5. determine the new bath Green matrix (self- 
consistency condition) 



g-\uj)^i:c{^) + G-\uj). 



(5) 



These steps are iterated to self-consistency. 

The Anderson model to be solved in step [5] is given by 



^And = ^clu+ 2J ^lm.,<y Aa^ma+zL ^'' [Aa'^'^'y + ^-C 



(6) 
where the operator a\^ creates an electron of spin a on 
bath-site /. The cluster Hamiltonian iJciu is obtained 
from the original Hamiltonian ([T]) by transforming to the 
reciprocal space of the super-lattice of clusters, and pro- 
jecting to the cluster. Writing the single-electron part of 
i?(k) as the matrix H(k), the single-electron part of _ffciu 
is given by 



Hc= / dkH(k), 



(7) 



The interaction terms are simply those of ([T]), restricted 
to the cluster. 



The Hamiltonian if (k) in the reciprocal space of the 
super-lattice {f } of clusters can be obtained by changing 
to the basis of operators 



-CDMFT 



(k)-E 



e '•'*^Cf+R,,^. 



(8) 



The resulting quantum cluster approximation is 
CDMFT. Alternatively, we can start from the operators 
in the reciprocal space of the lattice to obtain 



-DCA 



(k) = E' 



-ik(f+Ri 



Ri 



(9) 



Now we obtain the DCA. The choice of the operators 
in the two approaches differs just by local phase factors. 
In CDMFT this gaugei^ is chosen such that phases ap- 
pear only in matrix elements involving different clusters. 
Thus all matrix elements on the cluster are the same as 
in the original Hamiltonian. The price for retaining the 
original matrix elements on the cluster is a breaking of 
the translation-symmetry of the original lattice. DCA 
opts instead to retain this symmetry by distributing the 
phase change uniformly over the cluster-sites. The price 
for retaining translation-invariance is that the matrix ele- 
ments in the cluster Hamiltonian differ from those in the 
original Hamiltonian (coarse graining). In both cases, 
CDMFT and DCA, the eigenvalues of H(k) are identical 
to the eigenvalues of the non-interacting part of H 



III. HYBRIDIZATION SUM-RULES 

While the most general parametrization for the bath 
is given by expression Q r''^ we can always diagonalize 
the hopping matrix E among the bath-sites to obtain 



QaU{^u Vi}; Lo)=u + ^i-^c-Y,^' ^' 



UJ - £i 



(10) 



The hybridization matrix is then given by the tensor 
product of the vectors V/ , where 



Vu = E r«. 



rn (pi.m 



(11) 



and (pi are the eigenvectors of E with eigenvalues £i . 

To obtain sum- rules for the hybridizations, we write 
the inverse of the bath Green matrix as 



g-'iiv) = Sc(cj)+ / dk W + ^ - H(k) - Se(cJ 



Considering the limit a; ^ oo, expanding to order l/w^, 
using ([7]), and comparing to pU]) we find 



^V;V,^= f dkU^{k)- ( f dkn{k) 



(12) 



To illustrate this hybridization sum-rule we consider a 
representative set of examples. 



A. Single site 

We consider a d-dimensional lattice with hoppings i„ 
to the z„ n,*''-ncarest neighbors. For Nc = 1 we have 
H(k) = e^. Thus we find for the hybridizations 



E^.^ 



1 



{2ttY 7_ 



d^kei 






^n. 1^^ 



(13) 



where the integral is just the second moment of the den- 
sity of states, so that the last equation follows as in the 
recursion method J^ For a Bethe lattice of connectivity z 
with hopping matrix element t/^ the sum-rule reduces 



B. CDMFT 

We start by considering a linear chain with nearest 
neighbor hopping t and a three-site cluster Nc = 3. In 
the CDMFT gauge we have 



H(fe) 



1 e~ 


:iik 


1 


1 


g3.fc 1 






(14) 



so that He is the original single-electron Hamiltonian re- 
stricted to the cluster: 



H, = — 



3 

2^ 



7r/3 / 1 

dkU{k) = t \ 10 1 
'V3 V 1 



(15) 



The sum-rule p^ then is 



J2^iaA^ 




(16) 



i.e., only the sites on the surface of the cluster couple to 
the bath. 

The general CDMFT hybridization sum-rule (fT^ can 
be easily visualized: The integral over the Brillouin-zone 
of the cluster projects the single-electron part of the full 
Hamiltonian to the cluster (see eqn. ([7])). The matrix ele- 
ments of H^ are thus the two-step hoppings that are pos- 
sible on the cluster. Likewise the integral over the Hamil- 
tonian squared gives the second moments, only that here 
the intermediate site is not restricted to the cluster. Thus 
the sum-rule matrix is given by the second-order paths 
between cluster-sites that proceed via a site outside the 
cluster. This is illustrated in figure [TJ As a special case, 
for a single site we obtain the second equality in P^ . 

The vanishing of a matrix element in the sum-rule only 
implies that the corresponding matrix element of the 
bath Green matrix decays faster than l/uj for large uj. 
For a diagonal element, however, all terms in J^i yi,iyi,i 
are positive. Thus a vanishing sum means that all 
terms must be zero. Hence the sum-rule implies that 



a) 



b) 



c) 
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FIG. 1: CDMFT sum-rules for a one-dimensional 3-site clus- 
ter and nearest and next-nearest neighbor hoppings t and t" , 
respectively: a) E, |Vi,ip = t^ + t"\ b) E, Vi.iVi.a = «", 
and c) Ei l^'^i^ap ~ 2t"^ . The hybridizations are given by 
the two-step hopping processes that are lost when cutting the 
cluster out of the original lattice. 



cluster-sites that are so far in the interior that they can- 
not be reached by hopping from outside the cluster do 
not couple to bath and that all matrix elements of the 
bath Green function involving such a site i are given by 
G~j^{uj) = UJ + fj,— (H.c)ij for all uj. In that sense the bath 
hybridizes only to the surface of the cluster and we see 
that the hybridization-strength to these sites does not 
decrease for increasing cluster size Nc- 



C. DCA 

We start again by considering the 3-site cluster. In the 
DCA gauge we write 







H(fc) 



p Z /C p IK 



ik 



e^^ e 



e 



(17) 



Now He has translation symmetry, but the hopping ma- 
trix element is rescaled by sm{Tr / Nc) / {tt / Nc): 



He 



27r 



7r/3 o /o / 1 1 

d~kU{k)^-^^t 1 1 I . (18) 



-7r/3 



2tt 



1 1 



Since all matrices in ((T2l) are periodic, it is convenient 
to transform to fc-space. With V^k = X^i V'l^ic'-^^' /^/7^ 
and the coarse-graining factor r = 3\/3/27r we find 



/ , \Vl,K=±27r 



/3l 



2\,2 



{2-T/2-T^)t 



The hybridization sum-rule P^ is then, likewise, di- 



agonal in the cluster-momenta K 



(19) 



while all terms V/.k^/.k' mixing different cluster mo- 
menta vanish. As a special case, for a single site the 
above sum-rule is just the first equality in p3|) . Expand- 
ing £K+k around K, we find that for a d-dimcnsional 
system ^j |V;,kP decreases with cluster size as 1/Nc , 
while all cluster-sites couple with the same strength to 
the bathJ^ 



D. Discussion 

From the sum-rules we recover— that the individual hy- 
bridizations in CDMFT arc independent of cluster size, 

while for DCA they decrease with cluster size as A^c . 
Interestingly this means that for a d-dimensional system 
in CDMFT the overall coupling to the bath scales as 
Nc , while in DCA it scales a,s Nc ■ For non- 

local properties a DCA calculation is therefore expected 
to converge faster with cluster sizciii For a calculation 
where we represent the bath by discrete degrees of free- 
dom this decrease in hybridization strength does, how- 
ever, not help very much as we still need bath-sites to 
fit the hybridizations, even if they are small. With in- 
creasing DCA cluster size we thus have to parametrize 
Nc baths, one for each K. In CDMFT the situation is 
more fortunate, as the sum-rules imply that many hy- 
bridizations vanish and we only need to parametrize the 
coupling of surface-sites to the bath. 

The lack of translational invariance in CDMFT has two 
important practical implications. First, the full Green 
matrix has to be calculated, instead of just its diagonal. 
Second, when calculating local quantities, like the density 
per site, in CDMFT we have a choice of inequivalent 
sites, or we could consider the average over all sites. In 
a gapped system the best choice is the innermost site^^ 
however, in such a situation it might be better to do a 
straight Lanczos calculation with Nc + Ni, cluster sites, 
instead of using At bath sites f^^ 



IV. SYMMETRIES 

In the absence of spontaneous symmetry breaking the 
symmetries of the cluster (point-symmetries in CDMFT 
and additionally translation symmetry in DCA) arc re- 
flected in the Green matrix. In a symmetry broken 
state with long-range order, like an antiferromagnet or 
a charge-density wave, the symmetry of the Green ma- 
trix is accordingly lowered. To exploit this symmetry 
we introduce vectors on the cluster that transform ac- 
cording to its irreducible representations. We write these 
vectors as w/,y where / is the irreducible representation 



and V = 1 . . . Nj counts linear independent vectors trans- 
forming according to /. On an A^c-site cluster we can 
choose Nc such vectors that are orthonormal. Defin- 
ing the matrix W = i'^i^u) of these vectors, we can 
block-diagonalize the bath Green matrix: W^^^^W has 
blocks of dimension Nj corresponding to the different 
irreducible representations /. Since W^^^^^W is block 
diagonal for all cj, it follows from equation pTI)l . that 
W must also block-diagonalize the individual hybridiza- 
tion matrices V/ Vj . Therefore the hybridization vectors 
must transform according to an irreducible representa- 
tion: They can be written as V; = ^^ ^ij,!^ ^i,v for 
some irreducible representation /. If the V/ also had 
components wj ,y of a different irreducible representation 
J ^ I this would produce a hybridization matrix that 
could not be block-diagonalized. 

We thus find that the bath-sites can be arranged into 
sets corresponding to the different irreducible representa- 
tions. For fitting a block of the symmetrized bath Green 
matrix we need then only consider bath-sites of the re- 
spective irreducible representation. If the block is one- 
dimensional we can choose the corresponding hybridiza- 
tion real. An early example is the bonding-antibonding 
transformation introduced in reference |20|. 



A. CDMFT 

As an example we consider a linear cluster of 3 sites 
as shown in figure [51 The symmetry is C2 (see table U). 
Transforming to the basis vectors -wa,i = (|1) + |3))/V2 
and WA,2 = 1 2) of symmetry A (see table |l| and ws = 
(|1) — 13))/\/2, we find the transformed bath Green matrix 




WTg-iw = 



A bath-site of irreducible representation A contributes 
to the first block and has the same hybridization Va,i 
to the outer cluster-sites plus an independent hybridiza- 
tion parameter Va,2 to the central site. A bath-site of 
irreducible representation B contributes to the second 
block. For such a bath-site the hybridization to cluster- 
sites that arc related by mirror symmetry have opposite 
signs. Consequently, the hybridization to the central site 
vanishes. 

The situation is slightly more complicated when the 
symmetry group has irreducible representations of di- 
mension higher than one. The simplest example is 
the 2x2 cluster with C^y symmetry. With 'wai = 
(|1) + |2) + |3) + |4))/2, WB, = (|1) - |2) + |3) - |4))/2, 
and the pair w^.i = (|1) - |2) - |3) + |4))/2, we.2 = 
(|1) + |2) - |3) - |4))/2 we find that W^g-^W is diago- 



A •^v,,^^ 



B 



v„ 



^o-^# B 



FIG. 2: Hybridization of bath-sites of symmetry A and B to 
a 3-site cluster. As defined in table|IJ A is the unit representa- 
tion, so a bath-site of type A has the same hybridization V to 
all cluster-sites that are equivalent by symmetry. B is the an- 
tisymmetric representation, so the hybridization of a bath-site 
of type B to cluster-sites that are related by mirror symme- 
try have the opposite sign. Consequently the hybridization 
to the central site of a linear cluster with an odd number of 
sites vanishes in the B representation. 
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FIG. 3: Hybridization of bath-sites of symmetry Ai, B2, and 
i5 to a 2 X 2 cluster. For a given irreducible representation 
the absolute value of the hybridization to all cluster-sites is 
the same, while the signs are indicated in the figure. Non- 
trivial hybridizations corresponding to irreducible represen- 
tations A2 or Bi only appear for larger clusters. 
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\2))/V2 and y^e.2 = (|1) + |2) - 2|3))/V6. 



TABLE I: Character tables of the point groups Ci^,, C21,, C31,, 
and Cav 



nal with diagonal elements 



(wtg-iw)ii 


= Gil + 2^12^ + Gu 


(Wtg-lW)22 


= Gil - 25^2^ + G{3^ 


(Wtg-lW)33 


= Gii'-Gii 


(Wtg-lW)44 


= Gil - Gu 



A bath-site of symmetry Ai has the same hybridiza- 
tion to all cluster-sites while for a bath-site of symme- 
try i?2 the hybridizations have alternating signs: Vj = 
Vi(l,— 1,1,— 1). To realize the two-dimensional repre- 
sentation E we need two bath-sites li and I2 with de- 
generate energies ei^ = Si^ = ei and hybridizations: 
\l = Vi (1, -1, -1, 1) and Ml = Vi (1, 1, -1, -1). This 
is illustrated in figure [31 



B. DCA 

As an example for DCA we consider a 3-site cluster 
with periodic boundary conditions. The symmetry group 
is C3„ (translations and inversion). Hence we introduce 
the basis vector w^^ = (|1) -I- |2) -|- |3))/-\/3, correspond- 
ing to A: = 0, while the vectors formed by sin(27r/3) 
and cos(27r/3) give the E representation: -we,! = (|1) ^ 



w^g^w 




Gil G12 



In general bath-sites corresponding to the gamma point 
have the same hybridization to all cluster-sites, while 
those corresponding to fc = tt have alternating hybridiza- 
tions. For all other fc-points we need two degenerate 
bath-sites, with hybridizations VJ^.^ = Vi sin(fc/i) and 
Vi^^^ = Vi cos(fc^) to cluster-site fi. 



C. Cluster replica 

Instead of implementing the symmetry of the Green 
matrix as described above, one might construct the bath 
out of replica of the A^c-site clusterSfi^ For a two-site 
cluster this means that bath-sites come in pairs, with 
on-site energy e, hopping — t between the bath-sites, and 
hybridization Vn and V12 to the cluster as illustrated in 
figure m Diagonalizing such a bath-pair, we obtain one 
bath-site of symmetry A with on-site energy e^ = e — t 
and hybridization Va = {Vu + Vi2)/'\/2 and one bath-site 
of symmetry B with Eb = £ + t and Vb = (Vii — Vi2)/-\/2. 

To generalize this approach, let Hi, be the Hamilto- 
nian for an A^c-site cluster with general on-site energies 
and hoppings that respect the symmetry of the Green 
matrix. Furthermore let bath-site n hybridize to cluster- 
site i with Vni ■ These hybridizations are chosen symmet- 
ric under simultaneous symmetry transformations of the 
original cluster and the bath replica. Diagonalizing Hi, 
we obtain N^ eigenstates ipi with energy £;. These can 
be considered as bath-sites that hybridize to cluster-site 

i with Vl^i = J2n 'Pl.n ^n,i- 

Being the eigenstates of the Hamiltonian Hb, the ipi 
transform according to the irreducible representations of 



-Vii;:; 



FIG. 4: Replica of a 2-site cluster in the bath. Diagonalizing 
the Hamiltonian for the two bath-sites leads to a bath-site 
of type A (even representation) with energy ea ~ i — t and 
hybridization Va ~ {Vn + Vi2)/'\/2 and a site of type B with 
eB=e + tandVB = {Vn - Vi2)l\Pi. 



the symmetry of the Green matrix.— Therefore they can 
be written as hnear combination ipi = J^u oci;i ,v'w i .i' for 
some irreducible representation /. From this we can 
conclude that a cluster replica gives rise to TV/ bath- 
sites of symmetry /. By working with cluster replica 
we thus sacrifice the freedom of choosing the irreducible 
representations for the bath individually. Moreover, it 
is not straightforward to find a proper parametrization. 
Hi, must be chosen such that all accidental degeneracies 
can be lifted. For symmetries with higher-dimensional 
irreducible representations there will be, however, corre- 
sponding essential degeneracies. Moreover, working with 
cluster replica, we cannot fit the individual blocks of the 
bath Green matrix with the minimal set of symmetry- 
adapted parameters, but have to solve the optimization 
problem for the full bath Green matrix and all parame- 
ters. Thus, using cluster replica is less flexible than us- 
ing individual irreducible representations and it leads to 
a more complicated fitting procedure, in particular when 
considering large baths. 




FIG. 5: (Color online) Density as a function of chemical po- 
tential for a Hubbard model on the Bethe lattice with half 
band width D and U = 2 D in single-site DMFT with param- 
agnetic (para) and antiferromagnetic (AF) bath (blue squares 
and red circles respectively). The AF curve comprises a left 
and a right branch: the left one has been obtained by fix- 
ing the magnetization from 0.01 to 0.82 (from left to right) 
and determining the chemical potential. The right branch 
has been obtained by just decreasing the chemical potential 
from the half-filling value (p = U /2) down to the point where 
convergence is no longer found (/x = 0.6). The magnetization 
in the right branch is about 0.85. In the inset we compare 
the total energy of the AF and para phases as a function of 
doping. The AF phase is the stable one up to about 5 — 0.2. 
We have added S U/2 — EafJO) to E(5) to allow better com- 
parison with Fig. 1 of Ref. ^ and Fig. 8 of Ref. ^. 



ways metallic, the latter yielding a gap, which is over- 
estimated as we are using an antiferromagnetic bath to 
mimic the short-ranged correlations present in the Hub- 
bard chain. 



V. THE HUBBARD CHAIN 

We now discuss the one-dimensional Hubbard model, 
for which exact results are available from the Bethe 
ansata^^ and which has been studied with CDMFT us- 
ing Lanczosiiii^iii^ and QMC^^ as well as with the vari- 
ational cluster approximation.— Here going from single- 
site DMFT to a cluster description makes a qualitative 
difference: For a paramagnetic single-site calculation an- 
tiferromagnetism is completely suppressed, while on a 
cluster we will have short-ranged antiferromagnetic cor- 
relations, even if we impose a paramagnetic bath. 

This inclusion of antiferromagnetic correlations might 
well be the cause for the spectacular difference between 
the single-site and 2-site CDMFT calculations of the den- 

tas a function chemical potential reported in reference 
For illustration, in Fig. O we compare the density as 
a function of the chemical potential for a single-site calcu- 
lation with paramagnetic and antiferromagnetic bathj^ 
They represent two limiting cases, the former being al- 



A. Bath Green matrix 

For a one-dimensional lattice the CDMFT bath Green 
matrix simplifies drastically: When removing the A'^c-site 
cluster from the lattice, we arc left with two disconnected 
pieces. Thus a vanishing sum-rule p2p in one dimension 
means that the hybridization matrix vanishes for all fre- 
quencies. For the irreducible representations of the bath 
this means that bath-sites of symmetry A and B come in 
degenerate pairs: eai = £Bi and V^, = Vb; ■ Thus for the 
Hubbard chain with nearest neighbor hopping only the 
outer-most cluster sites hybridize with the bath and the 
bath parametrizations are identical by symmetryFi^ The 
evolution of the bath Green matrix clement Q^^ with 
cluster size is shown in Fig. [51 We find that the bath 
Green matrix elements hardly depend on cluster size and 
even for a single-site calculation the bath is already sim- 
ilar to that for a large cluster. 

In contrast, in the DCA we get a non-vanishing hy- 
bridization for each _R'-point of the cluster. This is shown 
in figurc[71 While the hybridization strength per /v-point 




((n-|j.)/r 

FIG. 6: (Color online) Plot of the CDMFT bath Green func- 
tion Q^-^ (uj) ~ {uj + fj,) on the real axis for linear clusters with 
nearest-neighbor hopping t. Increasing the size of the cluster 
the hybridization hardly changes. The plots show calculations 
for chains with U = 6t and n — 0.5f. 



decreases with cluster size, wc still have to parametrize 
all of them, possibly except for K ~ and tt, which 
almost vanish already for moderate cluster sizes. 



B. Fitting the bath Green matrix 

The most critical step in calculations with a finite bath 
is the determination of the parameters for the impurity 
Hamiltonian ^ . This is usually done by fitting the bath 
Green function on the imaginary axis: A fictitious tem- 
perature 1/(3 is introduced and the sum of the squared 
difference between Q and its parametrized version pop 
over the Matsubara frequencies up to some cutoff is 
minimizedi^ This procedure is fairly robust if the number 
of bath-sites is sufficiently large. For cluster calculations 
the number of effective bath-sites per fitted bath Green 
function can, however, be quite small. In such a situa- 
tion details of the fitting procedure are important and 
are accordingly discussed in the literature . ^^i^^'^'^i^^'^^i'^° 

To fit the Anderson parameters Vi and ei we use the 
distance function 
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where || • ||/ is the 1-norm for the block of irreducible 
representation / and N determines how strongly large 
Matsubara frequencies uJn are weighted. The distance 
function will only be finite in the limit Wcut — > oo if the 
sum-rules are fulfilled exactly, or if iV > 2. Thus for 
N < 2 increasing the cutoff emphasizes the sum-rules. 
The same is true for decreasing /?. If the distance func- 
tion is dominated by the large frequency asymptotics, the 
optimization mainly focuses on just the sum-rules. This 
means that in practice the Anderson parameters can be- 
come under-determined. Since the fictitious temperature 
is only used for fitting, while the calculations are actually 
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FIG. 7: (Color online) Plot of the DCA bath Green func- 
tions Q^^ {u}) — {u! + n — J dk e^-i-fe) '-'^ ^^^ ^^^^ ^^^^ ^'^^ linear 
clusters with nearest-neighbor hopping t. Increasing cluster- 
size {Nc =2, 4, 6, and 8 from top to bottom) the number 
of independent functions to fit increases. Note that K = 
and K = TT are associated with single bath-sites, while all the 
other A'-points need to be described by pairs of bath-sites. 
From the figures it is clear how the total spectral weight de- 
creases with increasing Nc- The plots show calculations for 
chains with U = 6t and fi — 0.5f. 
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FIG. 8: (Color online) Density for a 2-site cluster with [/ = 4i 
and /i = in CDMFT (left panel) and DCA (right panel) 
calculations with 12 bath-sites using different fictitious tem- 
peratures 1//3 but constant cutoff o^cut ~ 200t and A'' = 
and 1, starting from different Anderson parameters. For low 
fictitious temperature results are fairly independent of the 
starting points, while for larger temperature the Anderson 
parameters are essentially under-determined by the distance 
function and consequently the results of the supposedly self- 
consistent calculation strongly depend on the initial values. 



for T = there is a strong dependence of the physical 
quantities on the Anderson parameters. Hence for small 
/? the self-consistent results strongly differ for different 
initial Anderson parameters. This is illustrated in figure 
[51 Interestingly, the situation is opposite to calculation 
at finite temperature, where at higher temperature the 
physical quantities become less dependent on the details 
of the fitting^ 

To avoid an under-detcrmination of the Anderson pa- 
rameters, it is important to ensure that the features of 
the bath Green matrix close to the real axis are prop- 
erly weighted in the distance function. This was already 
pointed out in reference [l2- We find that a good com- 
promise between fitting the large- and small-frequency 
behavior is given by /? = 256/i, A'' = 1, and Wcut = 200 f, 
which, if not explicitly specified otherwise, is used in the 
calculations reported in this work. 



C. Convergence Avith number of bath-sites 

To check how many bath-sites we need to reach a satis- 
factory fit of the bath Green matrix, we consider a 2-site 
cluster for increasing A^^;,. As example we show in figure [5] 
the density for /x = 0. For CDMFT we find that we need 
at least 8 bath-sites to obtain a converged density. Also 
for DCA we obtain convergence for A'^f, = 8. Both these 
results translate to 4 bath-sites per non-vanishing cle- 
ment of the bath Green matrix. It is interesting to note 
that DCA converges to a density above the Bethe-ansatz 
result. This could be an artifact of the 2-site cluster, for 
which the coarse-grained hopping is larger than t, be- 
cause for periodic boundary conditions the hopping on 
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FIG. 9: (Color online) Density as a function of the number 
of bath-sites Nb for a 2-site cluster for U = At and /i = 0. To 
give a measure of the reproducibility of the results, we plot 
the densities of several converged CDMFT and DCA runs 
(multiple symbols at the same Nb). 
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FIG. 10: (Color online) Fraction of hybridization sum-rule 
for a 2-site cluster as a function of Nb- As more and more 
bath-sites are added, the exact hybridization sum-rule is ap- 
proached. Shown is the percentage of the sum-rule for irre- 
ducible representation A in CDMFT and K = Gm DCA. For 
irreducible representation B and A' = tt we find very similar 
behavior. 



the cluster and across the boundary add up. Averaging 
over different choices of boundary conditions or going to 
a larger cluster might improve the situationi^ 

Checking the hybridization sum-rule for the diagonal 
elements of the bath Green matrix, we find that the den- 
sity is already converged while the sum of the hybridiza- 
tions ^; I V; p is only at about 70% of its exact value. As 
shown in figure [10] we need to go to even larger baths 
to properly fulfill the sum-rule. We tried also to fit the 
bath Green matrix imposing the sum-rule, i.e. fixing one 
Anderson parameter per irreducible representation. We 
found, however, that for small A*";, this does not give par- 
ticularly good results, as it weights the large frequency 
behavior of the bath Green matrix too strongly, while for 
large A''f, it is not necessary. 

Going to larger clusters, we expect that we will need 



more bath-sites for a converged calculation. A notable 
exception is CDMFT for the linear chain with nearest- 
neighbor hopping only. As discussed above, in this case 
there are two identical baths coupling to one surface site 
each. Since the bath Green function that these baths 
have to fit is fairly independent of cluster size (see figure 
in]) we expect that the number of bath-sites needed for 
convergence is independent of cluster size. 



D. Convergence with number of cluster-sites 

We now analyze the convergence of CDMFT with the 
number of cluster sites N^. As before we focus on the 
density n, which is shown in figure [TT] for U = At and 
chemical potential fi ~ and /i = t. Considering the 
series of odd or even Nc separately, we see that with in- 
creasing cluster size the average density systematically 
approaches the exact result for the infinite chain. Inter- 
estingly, going from a cluster with an even number of 
sites Nc to A^c + 1 the average density hardly changes. 
For both chemical potentials already the smallest cluster 
gives a significant improvement over a single-site calcu- 
lation. For /i = the exact density is basically obtained 
for Nc = 2. For jjl = t convergence to the infinite-chain 
result is only reached at Nc — &- 

In figure [12] we show how the density versus chemi- 
cal potential curve for CDMFT calculations of increas- 
ing cluster size approach the exact result for the infinite 
Hubbard chain. We find that the closer we come to the 
mctal-insulator-transition the harder it gets to reach the 
infinite-size limit. This does not come as a complete sur- 
prise, as the self-energy is expected to become strongly 
/c-dependent at the Mott transition.— i^ 

We note that our results for Nc = 2 agree with figure 
4.7 of Ref. [ij- They are not compatible with figure 4.1 
of Ref. [H and figure 2 of Ref. [H. We have checked that 
our calculation is properly converged by starting from 
a number of different initial points, always converging to 
essentially the same density. To achieve this, for chemical 
potentials fi> 1.3 1 we increased f3 from 256/t to 512/t, in 
line with the trend shown in figure [5] We can, however, 
get significant variations in the density by putting re- 
strictions on the bath parametrization. Using, e.g., only 
six bath-sites, the CDMFT result happens to be closer 
to the one for the infinite system, in the vicinity of the 
Mott transitions^ This is shown in figure [T3l For small 
baths the calculated density is very sensitive on Nf, and 
can be either larger or smaller than the density for the 
infinite chain. For such small baths results will therefore 
critically depend on the fitting. We can, e.g., artificially 
"improve" the result by forcing a pair of bath energies 
to zero. Other restrictions on the bath-parameters in- 
stead move the densities further away from the Bethe 
curve. In all these cases we find that the restricted bath- 
parametrization results in a significantly deteriorated fit 
of the bath Green matrix. I.e., the sensitivity of these 
calculations to technical details merely shows the effects 
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FIG. 11: (Color online) Density for linear CDMFT clusters 
of increasing size Nc with f/ = 4i at /^ = (upper panel) and 
fj, = t (lower panel). All cluster calculations are for A^;, = 8. 
Circles denote the average density per cluster-site. (Green) 
Open squares are the individual cluster-site occupations. The 
size of the squares indicates how close the site is to the center 
of the cluster. To assess the reproducibility of the calcula- 
tion we show for each even value of Nc the result of at least 
two to seven runs with different choice of the starting set of 
Anderson parameters. The dotted line represents the exact 
Bethe-ansatz result. 



of an inadequate fitting of the bath Green matrix. By in- 
creasing Nb, the bath-parametrization improves and the 
calculated density converges, as shown in figure [T3] Nev- 
ertheless, differently from the behavior at smaller chemi- 
cal potential (cf . Fig. [9]) , the converged value is substan- 
tially smaller than the Bethe ansatz one. As this does 
not improve much with increasing Nc, we can conclude 
that the clusters are still too small to accurately cap- 
ture the behavior of the infinite system close to the Mott 
transition. 



VI. CONCLUSIONS 

The central problem of dynamical nrean-field calcula- 
tions with a solver that uses a small number of sites is to 
find a good parametrization of the bath. To address this 
problem for dynamical cluster approximations we have 
presented a systematic formulation for the bath degrees 
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FIG. 12: (Color online) Density n as a function of chemical 
potential n a,t U — At for linear CDMFT clusters of increasing 
size and Ni, — 8 compared to the exact result for an infinite 
Hubbard chain (Bethe Ansatz). The crosses give the results 
from figure 4.7 of reference Ha . The inset shows how the 
convergence of the density with increasing cluster size to the 
Bethe result becomes progressively slower close to the Mott 
transition, uba — ncDMFT jumps when the self-consistent so- 
lution changes sector {Ni,Ni). Close to these sector changes 
the results slightly depend on the initial conditions, i.e., there 
is a hysteresis between calculations increasing or decreasing 
fi. This is shown as multiple symbols for a given chemical 
potential. The upper plot shows the average density for 2-, 
4-, and 6-site clusters, the lower plot the density on the two 
central sites. 



of freedom. We have found sum-rules which allow to 
identify what hybridizations vanish and hence need not 
be parametrized at all. In addition the sum-rules can be 
used to check convergence for small baths. For the non- 
vanishing hybridization functions, we have introduced a 
bath-parametrization based on the irreducible represen- 
tations of the cluster Green matrix. In this approach the 
fitting of the bath sites is broken into independent fits 
of irreducible blocks of the Green matrix. This leads to 
a significant simplification of the fitting procedure which 
is particularly important when dealing with large baths. 
The symmetry-based approach should also benefit the 
variational cluster approximation (VGA) ^ where the de- 
termination of the parameters requires a Lanczos calcu- 
lation in each optimization step. 
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FIG. 13: (Color online) Average density as a function of the 
number of bath-sites A^;, for 2- and 4-site clusters with U = 4i, 
IJ, = 1.3 1, and (5 = 512/f. As for smaller chemical potential 
the density is converged for Nb = S (cf. figure 121) • We observe 
that for larger fj, unconverged calculations with small baths, 
e.g. Nb — 6, can give densities closer to the one of the Hubbard 
chain. 



As an application we have revisited the Hubbard chain. 
While this one-dimensional problem is the worst case sce- 
nario for DMFT, which is exact in infinite dimensions, it 
is technically the easiest case for CDMFT, because it 
requires only a minimal bath which is essentially inde- 
pendent of cluster size. This allowed us to study the re- 
sults of GDMFT using linear clusters of increasing size, 
extending previous work that was limited to 2- and 3- 
site clustersJ^ii^ Analyzing the density as a function of 
chemical potential, wc find that results significantly im- 
prove already going from a single-site DMFT to a 2-site 
cluster and become systematically better for larger and 
larger cluster sizes. Close to the Mott transition the con- 
vergence with Nc critically slows down, implying that 
the /c-dependence of the self-energy gets more and more 
important. 

The fortuitous independence of the GDMFT bath on 
the cluster size for the Hubbard chain is lost in higher di- 
mensions. Already for a cluster as small as 2 x 2 wc have 
to fit three functions, one of which belongs to a doubly de- 
generate E representation. In that case Nf, — 8 translates 
to only two effective bath-sites per bath Green function. 
For comparison, in our one-dimensional calculations we 
need at least 4 bath sites per bath Green function for a 
converged bath. Hence, using a Lanczos solver, it is vir- 
tually impossible to converge a 2 x 2 CDMFT calculation. 
This is illustrated in figure [TH 

The same is true for DCA, since the number of baths 
increases with cluster size, independently of the dimen- 
sionality of the cluster or the nature of the hopping. 
For zero-temperature cluster calculations it is therefore 
mandatory to move to impurity solvers that can handle 
large baths, e.g., DMRG.— For these calculations with 
large baths the efficient parametrization of bath and fit- 
ting of irreducible blocks will become even more impor- 
tant. 
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FIG. 14: (Color online) Fitting the bath Green function for 
a 2x2 cluster in CDMFT with U = 12t, n = 0, Nb ^ 8. The 
imaginary part of the block-diagonalized bath Green matrix 
(shown are only the three inequivalent blocks labeled by the 
irreducible representation) and its corresponding fit are plot- 
ted for two values of the exponent A'' (cf. section IV B|) . The 
quality of the fit is clearly very poor in both cases (for A'^=l 
compare full red/gray line with full black one and for A''=0 
compare dashed green/gray line with dashed black one). 



We would like to thank M. Capone for inspiring the 
present systematic treatment of CDMFT and for impor- 
tant discussions in the early stage of the work. M. Civ- 
elli kindly shared with us many details from his PhD 
thesis. A. Parola and S. Sorella provided us with the 
density curve for the Bethe Ansatz solution. We also ac- 
knowledge useful discussions with C. CastcUani. M. Hef- 
tier, M. JarrcU, S. Kancharla, G. Kotliar. A.I. Lichtcn- 
stein, I. Mazin, M. Potthoff, S. Sakai, and A. Toschi. 
S.G. thanks the Forschungszentrum Jiilich for hospital- 
ity. Calculations were performed in Jiilich on the JUMP 
computer under grant JIFF22. 



W. Metzner and D. VoUhardt, Phys. Rev. Lett. 62, 

324 (1989); M. Jarrell, tbtd. 69, 168 (1992); M. Rozen- 

berg, X.Y. Zhang, and G. Kothar, ibid. 69, 1236 (1992); 

A. Georges and W. Krauth, ibid. 69, 1240 (1992); M. Caf- 

farel and W. Krauth, ibid. 72, 1545 (1994). 

A. Georges, G. Kotliar, W. Krauth, and M.J. Rozenberg, 

Rev. Mod. Phys. 68, 13 (1996). 

M.H. Hettler, A.N. Tahvildar-Zadeh, and M. Jarrell, 

Phys. Rev. B 58, 7475 (1998); M.H. Hettler, M. Mukher- 

jee, M. Jarrell, and H.R. Krishnamurthy, Phys. Rev. B 61, 

12739 (2000). 

A.I. Lichtenstein and M.I. Katsnelson, Phys. Rev. B 62, 

R9283 (2000). 

G. Kotliar, S.Y. Savrasov, G. Palsson, and G. Biroli, 

Phys. Rev. Lett. 87, 186401 (2001); G. Biroli and 

G. Kothar, Phys. Rev. B 65, 155112 (2002). 

M. Potthoff, Eur. Phys. J. B 32, 429 (2003); M. Potthoff, 

Phys. Rev. Lett. 91, 206402 (2003). 

T. Maier, M. Jarrell, T. Pruschke, and M.H. Hettler, 

Rev. Mod. Phys. 77, 1027 (2005). 

A. Toschi, A.A. Katanin, and K. Held, Phys. Rev. B 75, 

045118 (2007). 

H. Hafermann, S. Brener, A.N. Rubtsov, M.I. Katsnelson, 

and A.I. Lichtenstein, JETP Lett. 86, 677 (2007). 

M. Jarrell, T. Maier, C. Hushcroft, and S. Moukouri, 

Phys. Rev. B 64, 195130 (2001); T. Maier, M. Jarrell, 

T. Schulthess, P. Kent, and J. White Phys. Rev. Lett. 95, 

237001 (2005). 

C.J. Bolech, S. S. Kancharla, and G. Kotliar, Phys. Rev. B 

67, 075110 (2003). 

M. Capone, M. Civelli, S. S. Kancharla, C. Castellani, and 

G. Kothar, Phys. Rev. B 69, 195105 (2004). 

M. CiveUi, PhD Thesis, Rutgers University (2006), 



arXiv:cond-mat/0710.2802 
" W. Kohn, Phys. Rev. 133, A171 (1964). 
^5 R. Haydock, V. Heine, and M.J. Kelly, J. Phys. C: Solid 

State Phys. 5, 2845 (1972). 
^^ Th. Maier, M. Jarrell, Th. Pruschke, and J. Keller, 

Eur. Phys. J. B 13, 613 (2000). 
^"^ K. Aryanpour, Th. Maier, and M. Jarrell, Phys. Rev. B 

71, 037101 (2005). 
^* G. Biroh and G. Kothar, Phys. Rev. B 71, 037102 (2005). 
^® E. Koch and S. Goedecker, Solid State Commun. 119, 105 

(2001). 
^" V.V. Mazurenko, A.I. Lichtenstein, M.I. Katsnelson, 

I. Dasgupta, T. Saha-Dasgupta, and V.I. Anisimov, 

Phys. Rev. B 66, 081104 (2002). 
^^ In case of accidental degeneracies, we can always sym- 
metrize the degenerate eigenstates. 

22 E.H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20 1445 (1968). 

23 B. Kyung, G. Kothar, and A. M.S. Tremblay, Phys. Rev. B 
73, 205106 (2006). 

2" M. Balzer, W. Hanke, and M. Potthoff, Phys. Rev. B 77, 

045133 (2008). 
2^^ G. Sangiovanni ei al, Phys. Rev. B 73, 205121 (2006); 

G. Sangiovanni et al, Phys. Rev. Lett. 97, 046404 (2006). 
2"^ P.G.J, van Dongen, Phys. Rev. B 54, 1584 (1996). 
2^^ R. Zitzler, T. Pruschke, and R. Bulla, E. Phys. J. B 27, 

473 (2002). 
2** B. Kyung, S.S. Kancharla, D. Senechal, A.-M. S. Tremblay, 

M. CiveUi, and G. Kotliar, Phys. Rev. B 73, 165114 (2006). 
2^ M. CiveUi, M. Capone, A. Georges, K. Haule, O. Parcol- 

let, T.D. Stanescu, and G. Kotliar, Phys. Rev. Lett. 100, 

046402 (2007). 
3° Y. Z. Zhang, and M. Imada, Phys. Rev. B 76, 045108 

(2007). 



31 
32 



12 

A. Toschi, private communication. Phys. Rev. Lett. 93, 246403 (2004); S. Nishimoto, 

M. Jarrell, private communication. F. Gebhard, and E. Jeckelmann, J. Phys.: Condens. Mat- 

M. Civelli, M. Capone, S.S. Kancharla, O. ParcoUet, and ter 16, 7063 (2004); M. Karski, C. Raas, and G. Uhrig, 

G. Kotliar, Phys. Rev. Lett. 95, 106402 (2005). Phys. Rev. B 72, 113110 (2005); D. Garcia, E. Miranda, 

T.D. Stanescu, M. Civelh, K. Haule, and G. Kothar, An- K. Hallberg, and M. Rozenberg, Phys. Rev. B 75, 121102 

nals of Physics, 321, 1682 (2006). (2007). 
D. Garcia, K. Hallberg, and M. Rozenberg, 



